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Abstract 

We have estimated contributions from the lossy and inductive vacuum chamber components 
to the broadband impedance of the NESTOR storage ring by using analytical formulas. As was 
expected considering the small ring circumference (15.44m), the main contributions both to the 
longitudinal impedance Zn/n and the loss factor ki oss come from the RF-cavity. 

Cavity impedance was also estimated with CST Microwave Studio (CST Studio Suite™ 2006) 
by simulating coaxial wire method commonly used for impedance measurements. Both estimates 
agree well. The upper limit of impedance of elliptic holes in the vacuum chamber of dipole magnet 
were also obtained with this approach. 

We have also evaluated the bunch length in NESTOR taking the conservative estimate of 3 Ohm 
for the ring broadband impedance and have found that the bunch length a z = 0.5 cm could be 
obtained for the designed bunch current of 10 mA and RF-voltage of 250 kV. 
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Introduction 

At the design stage of any storage ring it is customary to study the effects of beam 
interaction with vacuum chamber because these effects determine, for the most part, the 
bunch parameters and limit the stored beam current. Since pioneer work by Sessler and 
Vaccaro [1] the effects of beam interaction with a vacuum chamber (or a beam-pipe) are 
treated in the frequency domain in terms of coupling impedances: longitudinal impedance, 
Z\\(u) and transverse impedance Z±(u>). Generally, the beam current is characterized with 
higher moments which describe the dynamics of charge distribution in the bunch, and one 
has to match to all moments m=l,2,3... the corresponding impedances Z™(oj) and Z™(u) 
(see, for example, [2]). In practice only impedances with m=0,l which describe rigid bunch 
movement are studied in order to optimize the design of beam-pipe components and to con- 
sider the forthcoming beam instabilities and their possible cures. In this paper we consider 
the longitudinal broadband impedance which determines the bunch length in the NESTOR 
ring. 

The longitudinal impedance is related to longitudinal beam component (m=0), which for 
a point charge is I z = Iq8{x — xi)5(y — y\) exp(— i^z), and can be expressed in the following 
way: 

oo 

Z„(w) = ZT=\uj) = -i / E a e- iue '*dz , (1) 

— oo 

where E z is the longitudinal component of the electric field along the beam axis, c is velocity 
of light. Minus in Eq.Q imply that the induced field is retarding (shifted in phase by 180°). 

For convenience, one usually consider broadband (BB) ring impedance, that corresponds 
to all low-Q components of the beam chamber, and narrow-band resonances originated from 
high-Q elements, first of all, RF-cavities. RF-cavities along with fundamental (accelerating) 
mode TMoio have a whole spectrum of higher order modes (HOM's) which can drive coupled- 
bunch instabilities (CBI) [3]. Cavity HOM's as a possible origin of CBI in the NESTOR 
were discussed in the previous paper jl]. 

Longitudinal broadband impedance is usually normalized to a harmonic number n = 
uj/uo, where uo is a rotation frequency, and is expressed as Z\\/n. For many beam-pipe 
elements Z\\{oj) is inductive, so it is supposed that Z\\(cu)/n is constant in a wide frequency 
range. Zu/n is used for estimation of single-bunch instability thresholds with help of half- 
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empirical criteria obtained in low-frequency approximation (u> <C c/a z , where a z is a bunch 
length). Obtained in this approach, Keil-Schnell-Boussard criterion [5] relates Zn/n to the 
microwave instability threshold: 



n 



27ra(E /c)(5| 

Ipeak 



< -"^r'~' E . ( 2 ) 



where a is the momentum compaction factor, Eq/c is the electron beam energy and Se 
is the relative beam energy spread. I pea k represents the bunch peak current which for the 
gaussian bunch of length a z is I pea k = I av y/27rR/a z: where I av is the average bunch current 
and R is the average radius of the ring. This approximation is good for a medium bunch 
length: o z < b, where b is the characteristic transverse size of the beam pipe. For short 
bunches (a z <C c/co r ) one has to use in Eq.Q instead of the low-frequency impedance Z\\/n 
the effective impedance (Z\\/n) e ff, which is averaged over bunch frequency spectrum, thus 
taking into account that at high frequencies the bunch does not interacts with the low- 
frequency part of broadband impedance. For short gaussian bunches a simple relation can 
be used: (Z\\/n) e ff ~ 2(co r a z /c) 2 ■ \Z\\/n\ 

The real part of BB impedance determines beam energy losses. Parasitic energy loss 
per turn due to bunch-environment interaction is proportional to the bunch charge squared 
(AE = —ki oss q 2 ), where factor ki oss is a loss parameter which for gaussian bunch can be 
presented as follows: 



hoss = 7r- / ReZ»(u)exp( ^)dz , (3) 

Z7T / C z 



Below the cut-off frequency of a beam-pipe, u cut - ff, broadband impedance is an additive 
value and the sum of contributions from various vacuum chamber components (so called 
"impedance budget") is usually evaluated in order to minimize this value. For the frequency- 
dependent components (like resistive wall) impedance is evaluated at the roll-off frequency, 
^roil-off (the frequency at which the bunch spectrum power density |p(w)| 2 is one half of its 
peak value). For gaussian bunch ui ro ii-off = cy/2/a z . 

Total parasitic energy losses have to be considered in designing of heat-removing circuits 
and in developing of the RF-system as a whole. The latter is especially important for low- 
energy rings like NESTOR, in which the parasitic losses exceed synchrotron radiation losses 
and define the synchronous phase of the beam. 

In this paper we have evaluated contributions from various beam-pipe components to the 
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ring longitudinal BB impedance. Because of unavailability of 3D time-domain codes like 
GdfidL[6j or new version of CST Studio Suite™ [7j, which can directly calculate Z\\/n, we 
had to restrict ourselves to analytical approaches. In most cases the latter give satisfactory 
results. We have also used the available version CST Studio Suite™ 2006 to evaluate the 
impedance of the elliptical holes in vacuum chambers of dipole magnets, which couldn't be 
estimated analytically. 



1. The BB impedance of the NESTOR. Analytical approach. 

The next ring components can give a substantial contribution to the ring broadband 
impedance: 

- resistive wall, 

- RF-cavity; 

- pumping slots in dipole vacuum chambers; 

- beam pick-up's; 

- holes in the laser-electron beam crossing chamber; 

- RF-liners (bellows) 

- welding joints. 



1.1. The resistive wall impedance 

The resistive wall coupling impedance of a beam-pipe of elliptical cross section with major 
and minor semi-axes a and b can be evaluated with the following equation [8]: 

l+j)G (q), (4) 



n 2b 2-nR 

where L is the pipe length, 5 = (2cp/uZ Y^ 2 is the skin depth of the wall material whose 
specific resistance is p, Z is the impedance of free space. Go(q) is a function of the parameter 
q — (a — b)/(a + b) and is calculated via elliptic integrals and Jacobi elliptic functions. For 
NESTOR g=0.49, function Go(q) ~ 1 and the resistive wall impedance is given by a simple 
equation: Z^ w = 1.35(1 + i)^/n [Ohm]. For the bunch length a z =l cm (u) ro ii-off ~ 25GHz) 
the resistive wall impedance ZP W /n amounts approximately to 0.13 Ohm. 
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The estimate of kf^ according to Eq.([3]) gives 0.06 V/pC for bunch length a z =l cm. 



1.2. The RF-cavity 

The contribution from the RF-cavity to the ring broadband impedance can be estimated 
by using the following equation [H]: 



where Ui and (R/Q)i are the resonance frequency and R/Q factor for the cavity mode with 
mode index i. The sum is taken over all HOM's with < u cut - ff- HOM's parameters were 
calculated with ANSYS code [TO] using approach described elsewhere [I]. The estimation 
gives Zf OM /n=lA0 Ohm. 

Cavity HOM's give a significant contribution to the total loss factor. This contribution 
was estimated as sum of energy losses over all cavity HOM's: 



It gives kJ^° M =0A8 V/pC for bunch length a z =l cm and 0.54 V/pC for a z =0.5 cm. 
Energy loss at fundamental mode doesn't depend on bunch length in this range of a z and 
amounts to 0.54 V/pC. 

1.3. Beam position monitors 

The beam position monitor (pick-up) represents four electrostatic electrodes (buttons) 
placed on the inside surface of the beam-pipe and separated from the latter by narrow 
annular slots. The beam-pipe cross section at pick-up location is given in Fig. [T] Impedance 
of one button can be estimated as impedance of an annular slot in the elliptic beam-pipe. 
The longitudinal impedance of of a hole in an elliptic beam-pipe in static approximation 
(hole dimensions are small compared with the wavelength) is given by [TT] : 




(5) 





Z\ 



hole 



juZ 
8vr 2 c ' 
2K(k) 



i tyv - x) ■ Ql(v) 



(7a) 
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2 + b 2 sinh 2 m + sin 2 v 



Qo(v) 



k' 



(7b) 
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FIG. 1: Beam pipe cross section at pick-up location 



where ip v and \ are the susceptibility and polarizability of the hole; v is an azimuthal elliptic 



coordinate; u is related to ellipse semi-axes via: coshw = a/\/a 2 — b 2 ; dn(v,q) is Jacobi 
elliptic function of the argument v = 2K(k)v/iv)] K(k) is the complete elliptic integral of 
the first kind; k' = yl — k 2 and q is defined above. The azimuthal elliptic coordinate of the 
button center is v — 75.3° and Qo(v)=2.58. 

The susceptibility ip and polarizability x of a narrow annular slot (w = r ext — r int r ext , 
where r ext and r int are external and internal radii of the slot) in a thin wall are given by 
H2I: 



X 



2 2 
n r ext r int 



\n(32r ext /w) - 2 

--lX 2 W 2 (r ext + r int ) 



(8a) 
(8b) 



Calculations give Zf PM jn — 1.1 • 10~ 2 Ohm for one pick-up (4 buttons). It should be 
noted that the value of pick-up impedance obtained as a difference of impedances of two 
circular holes with radii r ext and ri nt |T3j is four times less than the value given above. 

For estimation of the real part of the pick-up impedance we used the formula obtained 
for a circular beam-pipe [H] and modified in accordance with Eq.([7]): 

Hole .{"V + X 2 ) ■ QI(V) 



ReZ, 



96vr 3 [b 2 + (a 2 -b 2 )sin 2 v] 



(9) 
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where Qo(v) is defined by Eq.Q. The loss parameter is given by: 

(10) 



.hole = ZoV^ (jg + X 2 ) • Q&V) 

vloss 167r%5 b 2 + (a 2 - b 2 ) sin 2 v 



Calculation gives k^ s M = 6.8 -10 5 • a z 5 [cm] V/pC for one pick-up. 



1.4. Welding joints 

The adjacent sections of the stainless-still beam-pipe in NESTOR are joined together by 
welding for which operation there are membranes at section ends. After welding between two 
joined sections an annular gap ~ 1.5mm wide and ~ 3mm deep is formed. The contribution 
of this gap to the ring broadband impedance has been estimated with the expression |15j : 

Z \\ aP = upZpgA 
n 2nca e ff 

where g and A are the gap length and depth, correspondingly, a e ff = \fab~ is the effective 
radius. We estimated Z?, ap /n to be ~ 4.8 • 10 -3 Ohm per joint. 



1.5 Bellows liner 



To match the environment seen by the beam in bellows with circular cross section to 
that of the elliptic beam pipe the RF-liners will be used. The impedance of such a liner 
can be estimated as impedance of many identical slots which are equally spaced around the 
azimuthal perimeter [Xlj : 

yliner • ry ( i \ 

z || ]Z r]{i) v - x) 



n (a + b) 



G (u ) (12a) 



47r J [stnh 2 u + sin l vy-i l 
o 

where rj = P/2ttRC, P is the number of slots, C is the circumference of the beam pipe cross 
section. For 22 slots, 2mm wide and 80mm long, the calculations give Z l ™ er /n = 2.1 • 10~ 4 
Ohm for one liner. 
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2. The broadband impedance of the NESTOR ring. Simulations. 



For some beam pipe elements the BB impedance couldn't be estimated analytically. 
Such are pumping slots in the vacuum chamber of dipole magnet, which have transverse 
dimensions comparable to those of the beam pipe. We have tried to estimate their impedance 
by simulating coaxial wire method, widely used for impedance bench measurements [16J, 
with available code pack CST Studio Suite™ 2006. This version includes transient solver 
that can be used for S-parameter calculation. 

The basic concept of this method relies on substituting the beam by a thin wire and thus 
simulating the fields of ultrarelativistic beam on the beam pipe wall by the propagation of 
TEM mode in the transmission line so formed. In the transmission line framework a single, 
lumped wall impedance Zy? can be expressed in terms of S-matrix coefficients in a following 
way: 

Su = or Z +7 and S21 = op 2 f C 7 ' ( 13 ) 

zix c -+- Zi\y ^tt c -+- Zj\y 

where R c is characteristic impedance of the transmission line. In principle either coefficient 
can be used for impedance determination, but the transmission coefficient is applicable in 
more general configurations and is usually preferred in bench measurements. 

The question to what degree the simulated impedance Rw approximate the BB 
impedance is raised from time to time [TTJUS]- A good agreement between calculations and 
measurements has been shown for a simple cavity-like structure with a thin wire (d=0.75 
mm) [TT] . 

So as the first step we applied this approach to estimate the broadband impedance of the 
RF-cavity and to compare it to that obtained with analytical approach. 



2.1 The RF-cavity 

The simplified cavity model (without ports) used in calculations is shown in Fig. [2] The 
calculated S- parameters in the frequency range 1 -j- 5 GHz presented in Fig. [3] 

The analysis of the electromagnetic field distribution at resonance frequencies reveals: 
1. The field patterns at resonance frequencies (minima in S21) i n peripheral region are 
similar to that of cavity modes. In particular, first resonance corresponds to TMqio mode 
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FIG. 3: Scattering matrix parameters: a) Sn (reflection) and b) S21 (transmission) 

but the mode frequency is shifted to higher frequencies. The field in the output beam pipe 
is zero i.e. the cavity acts clS db notch filter. 

2. At frequencies corresponding to minima in Su the electromagnetic wave doesn't excite 
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the cavity (field in peripheral region is zero) and the structure behaves as an usual coaxial 
transmission line. 



Cavity impedance obtained from with Eq.(13) is presented in Fig. 4. 
14-i 




2 3 
Frequency,GHz 

FIG. 4: Cavity impedance. 



Contribution to the ring broadband impedance from cavity HOM's was calculated with 
the following expression: 



7HOM 



/n 



1 



Z w (oo)—duj 



id 



(14) 



Bottom limit of integration u c = 2irf c was taken in the minimum between the first and the 
second resonances (it is shown by an arrow in Fig. Ill) . 



The obtained value Z¥ /n=1.5 Ohm well agrees with the estimate obtained with ana- 
lytical approach (see 3.1). It doesn't include contribution from propagating modes because 
^cut-off— 5.9 GHz for the cavity elliptic beam pipes. 
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2.2 The dipole chamber 

The cross section of the main vacuum chamber in the dipole magnet is similar to that 
presented in Fig. 1. The structure of the chamber is presented in Fig. [5j The main chamber 
is pumped with an ion pump, mounted on antechamber, through 16 pumping slots of elliptic 
cross section (22x14 mm). Two chambers of total number of four have additional holes for 
electron beam injection or X-ray beam extraction. The simulated model, which presents the 
"straightened" segment of the chamber with a single hole, is presented in Fig. [6] 




FIG. 5: A vacuum chamber in the dipole magnet (bottom half) 




FIG. 6: The model of the chamber fragment with a single hole used in simulations. 

To control the accuracy of obtained results we have made parallel simulations for similar 
elliptic segment without a hole. Impedance of the last is zero. Simulations have shown that 
obtained values of Sn are less than 10 -4 and are comparable to the level of numerical noise. 
The number of mesh points N c varied in the range 40000-^400000, at higher values of N c 
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the simulations became unstable. Calculations were performed on the typical office PC and 
were CPU-bound, so we had no opportunity to improve simulated model. 

Nevertheless, we have evaluated the upper limit of hole impedance by using strong az- 
imuthal dependance predicted by Eqs. ([7| for impedance of a small hole in the elliptic beam 
pipe with a substantial eccentricity. The elliptic hole (22x14 mm) in the elliptic pipe (79x27 
mm) was placed at various azimuthal angles perpendicular to the tangent plane. The thick- 
ness of the wall in the hole center was constant. The broadband impedance was calculated 



with Eq.(14), averaging was done over the range of 4- uj max . The obtained dependence of 
Zu ole /n on elliptic azimuthal angle v is presented in Fig. 7 
One can see that for azimuthal angles v < 30° the calculated impedance is nearly con- 




0,000 L-i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— ' 

10 20 30 40 50 60 70 80 90 



Azimuthal angle v, deg 

FIG. 7: The hole impedance versus a hole azimuthal angle. 

stant. It can be explained by the fact that at low angles the hole impedance decreases to 
values, which are below the level of numerical noise. The figure also shows that the single 
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hole impedance Z^° le /n >0.003 Ohm can be quite reliably deduced from S'n-parameter in 
coaxial wire simulations. For reasons given above the value 0.05 Ohm can be considered as 
the upper limit of the impedance of the dipole chamber (16 holes). One should also keep in 
mind that such an additive estimate is valid only in low-frequency approximation. 

3. The BB impedance of the NESTOR ring. Results. 

The obtained results are summarized in Tables 1,2. As follows from the tables for bunch 
length considered (<r z =0.5, 1.0 cm) the main contributions both to longitudinal broadband 
impedance and loss factor come from the RF-cavity. From a comparison of the loss factor 
data for two bunch lengthes one can see how the contributions from other components are 
increased with bunch shortening while the contribution from the RF-cavity isn't changed 
essentially. 



Table 1. Longitudinal BB impedance budget 



Component 


Number 


L, nH (1 unit) 


\Z\\/n\,n (total) 


RF-cavity 


1 




1.40 


Resistive wall 


1 


1.07 


0.13 


Dipole chamber 


4 


< 0.4 


< 0.2 


BPM 


4 


0.09 


0.044 


Welding joints 


8 


0.04 


0.04 


RF-liners 


4 


0.002 


0.0009 


Total 






< 1.65 



Table 2. Loss factors of beam pipe components 



Component 


Number 


hot,y/pc 


a z =lcm 


a z — 0.5cm 


RF-cavity 


1 


1.02 




1.08 


Resistive wall 


1 


0.06 




0.17 


BPM 


4 


0.0003 




0.0088 


Total 




1.08 


1.26 
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It should be mentioned that for the time being we have no estimates for two components: 
i) the beam-pipe section for the crossing point of the electron and laser beams; ii) the 
injection section (strip-line inflector). The first element at the commissioning stage will be 
replaced with a straight section. The injection section is essentially non-symmetric, it gives, 
presumably, a substantial contribution to broadband impedance and so it requires to be 
studied with 3D time-domain codes. 



4. Turbulent bunch lengthening in the NESTOR storage ring 

The turbulent bunch lengthening is observed at bunch currents higher some threshold 
value and is explained by coherent mode coupling arising through overlapping of mode 
frequencies at high bunch currents [19J . This effect displays itself as a longitudinal instability 
(microwave instability) and is accompanied by increase of the beam energy spread. 

The threshold bunch current 1^ for a given value of Z\\/n can be obtained with Eq.(|2]). 
Above threshold the bunch length can be estimated with the equation [20J: 



n 



/iKsin$< 



1/3 



(15) 



where h is the harmonic number, V c is the amplitude of RF- voltage and $ s is the synchronous 



phase of the beam. One can see from Eq.(15) that the bunch length doesn't depend on the 



electron beam energy above the instability threshold. 

We estimated bunch lengths for a reasonable range of RF- voltages. NESTOR parameters 
used in these calculations are presented in Table 3, the obtained results are given in Table 
4. 

In our calculations we have taken Z» / n=3 Ohm, that is almost two times the value given in 
the table, in order to make compensation for the contributions from beam-pipe components 
not considered in the paper, first of all, the inflector. The values of the energy spread 5^ s 
and bunch length af s presented in Table 4 were calculated with DECA code [21], which 
takes into account the effect of intrabeam scattering (IBS) - electron-electron scattering 
in a bunch. The threshold currents obtained with Eq.Q for these IBS-corrected bunch 
parameters are below the goal value of 10 mA per bunch, so the bunch will be lengthening 
in the ring. For the designed RF-voltage of 250kV estimation gives af l =0.5 cm. 
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Table 3. NESTOR ring parameters 



Parameter 


Units 


value 


Average radius, R 


m 


2.46 


Harmonic number, h 




36 


Compaction factor, a 




0.01 


Bunch current (maximal), I av 


mA 


10 


Ring broadband impedance, Z\\/n 


Q 


3.0 


Synchronous phase, $ s 


deg 


89 



Table 4. Bunch length in NESTOR 



V c , kV 


E =60 MeV 


£ =250 MeV 


erf 1 , cm 


8f 


af s , cm 


III mA 


8f 


af s , cm 


III mA 


10 


0.17 


10" 


-2 




1.32 


7.7 


0.47 


10" 


-3 




0.71 


1.2 


1.46 


50 


0.19 


10" 


-2 




0.67 


5.1 


0.53 


10" 


-3 




0.36 


0.8 


0.85 


180 


0.21 


10" 


-2 




0.40 


3.7 


0.59 


10" 


-3 




0.21 


0.6 


0.56 


250 


0.22 


10" 


-2 




0.34 


3.4 


0.60 


10" 


-3 




0.18 


0.5 


0.50 



At the commissioning stage of the NESTOR project the RF-amplifier with output power 
of 1 kW is available. It provides the accelerating voltage of 50 kV and the stored beam 
currents of up to 100 mA. We have no equipment to control a ring filling pattern, so all 36 
RF-buckets will be filled up, and a maximal bunch current will be decreased down to 3 mA. 
For this case estimation gives af l —0.6 cm. 

5. Conclusion 

The conservative estimate of the longitudinal broadband impedance of the NESTOR ring 
(Zn/n—3 Ohm) confirms the possibility of obtaining the bunch length of 0.5 cm with the 
goal bunch current of 10 mA for the designed RF-voltage of 250 kW. At the commissioning 
stage we can obtain 3mA in a 0.6 cm bunch with the RF-voltage of 50 kV. 
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